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ABSTRACT 


The multi-loop Lead Angle servo loop is utilized as a analog computing 
loop that mechanizes the range-bearing torpedo fire control equation. 
This thesis applies several new servo analvsis techniques, developes 
system component transfer functions, analyzes individual loops, and 
examines qualitative stability and response of the overall system. 

No response or error criteria was specified. Based on the author's 
experience and advice from submariners a minimum error criteria is 
derived. This thesis preposes several modifications of the origional 
system in order to improve performance, points out several difficiencies, 


and recommends areas of further study. 
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Definition 


When used alone, B denotes 

the relative bearing to the 
present position of the target 
and is measured in the Master 
own sh{p horizon system. 


The difference between 
Relative Target Bearing and 
Relative Tube Bearirg, mea- 
sured in the horizontal plane 
of the own ship horizon system. 


The angle formed by the 
intersection of the line from 
own ship's present position 

to target's present position 
with the line from torpedo's 
pseudo or real position to 
target's future position. 
Positive angles are measured 
clockwise from the line connect- 
ing own ship and target, i.e., 
the present relative bearing 
line. Prior to event lh, and 

for all Preset problems, 

point lh is the pseudo position. 
For Wire Guided problems, 

after event h, the dead reckoned 
torpedo position is used, 


Relative bearing of the torpedo 
tube measured in an own ship 
horizon system. 


Angle between the horizontal 
projection of the missile or 
torpedo speed vector and the 
present bearing line from own 
ship to the missile or torpedo, 
measured in the horizontal 
plane of a missile or torpedo 
horizon system. 
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Computed Correction 
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Definition 


Angle between the target speed 
vector am the present besrirg 
line frcm own ship to target, 
measured in the target horizon 
system, 


Course of own ship measured 
in. the Master system. 


Course of the torpedo or 
missile measured in the 
horizontal plane of a north 
oriented horizon system at 
The torpedo or missile. 


The course of the torpedo 
measured in the horizontal 
plane of a north oriented 
horizon system at point }). 


The course of the torpedo 
measured in the horizontal 
plane of an own ship oriented 
horizon system at point );. 


Course of the target measured 
in the target north oriented 
horizon system. 


The time rate of change of 
true target bearing. 


The rate of change of own 
ship course with respect 
to real time. 


The rate of change of target 
course with respect to time, 


The time rate of change of 
the distance traversed by 
own ship, measured in the 
horizontal plane of the 
Master own ship horizon 
system. 





Svmbol 


DMht 


Fe 


Fo 


e(G) 


sq(G) 


ce (Hm) 


ec(Hm) 
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Name 


Target Speed 


Frequency, cycles per 
second 


Closed loop transfer 
function 


Open loop transfer 
function 


Gyro Angle 


Error in Gyro Angle 


Angular Spread 


Gyro Angle at Event 6 


Length of Torpedo Path 


Computed Length of 
Torpedo Path 


Error in Computed Length 
of Torpedo Path 


Run to Burst Command 


Correction to Torpedo 
Path Length 
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Definition 


The time rate of change of 

the distance traversed by 

the target, measured in the 
horizontal plane of the target 
horizon system. 


Angle between vertical 
plane thru line of fire and 
vertical plane thru desired 
missile centerline measured 
in horizontal plane. (When 
used alone only.) 


A manual correction to gyro 
angle to initiate an angular 
spread of torpedo paths. 


The torpedo gyro angle at 
event 6. 


Actual distance traveled by 
torpedo from own ship to 
target. 


Manual correction to torpedo 
path length. 


A bias added to torpedo path 
Length to compensate for 
torpedo acceleration to 
required speed. 














Symbol Name Definition 
sq (Hm) Enabling Run Offset Correction used to shift 
enabling point along the 
torpedo track. 
Hm26 Total Run to Burst The total theoretical length 
of the torpedo path from event 
2 to event 6. 
(Obsolete symbol Hmt) 
Hm6 Distance from Torpedo to The straight line distance 
Point 6 from the torpedo's present 
position to point 6. 
Hm7 Remaining Run to Burst The total distance the torpedo 
must run from its present .....- 
J position to its bursti-n¥ position 
(point 7). | 
Hvm Torpedo Run Depth Vertical distance the torpedo 
is set to run below the surface 
of the water. 
Hvymeg Depth Difference Vertical distance between the 
torpedo tube and the missile 
or torpedo. (Hvm - Hvg) 
H97 Distance from event 
2 to event 7 
H67 Linear Spread Along Tar- The distance along the target track 
get Track between the future target position : 
(point 6) and the advance point (point7). 
4 
Ly2 Latitude of Firing The latitude of own ship at the time 
of fire. 
Lyo Latitude of Own Ship 
Lyp Latitude of Proofing The latitude at which the torpedo : 
was proofed, | 
. 
Lyt Latitude of Target 
k One Thousand ohms | 
K Gain 





Symbol 


Pdn 


Pdo 


PyoTo- ee 


J (R) 


c(Rh) 


J (Rh) 


Rh2 


Name 





Own Ship Horizontal 
Motion 


Athwartship Parallax 
Displacement 


Centerline Parallax 
Displacement 


Reots of equation 


Present Slant Range 


Computed Correction to 
Range 


Horizontal Range 


Computed Horizontal Renge 


Horizontal Range Correction 


Horizontal Range to Target 
at Event 2 


Definition 


The horizontal distance made good 
by own ship during a given inter- 
val of tine. 


The component, in the deck re- 
ference plane and perpendicular 

to own ship centerline, of the total 
distance from the origin of the 
Master system to the origin of the 
weapon system. (The center of the 


torpedo gyro or the missile platform. ' 


Pdn is measured in the Master own 
ship deck system. 


The corponent, along own ship center- 
line, of the total distance from the 
origin of the Master system to the 
origin of the weapon system. (The 
center of the torpedo gyro or the 
missile platform.) Pdo is measured 
in the Master own ship deck system, 


R, when used alone, is the 
present distance between 
own ship and target, mea- 
sured along a straight line 
in the Master system. 


The present horizontal range 
from own ship to the tarret's 
present position as seen by 
radar or as computed from 
sonar data, measured in the 
Master system, 


Computed correction to 
horizontal range. 


The range, measured in 

the Master system, from 
own ship to target at time 
of fire (event 2). 





T26 


Um 


j (Um) 


Nare Definition 
Horizontal Range to “each Tre range, measured in 
Point the Master system, from 


own ship present cosition 
to the torpedo reach point, 
potat 3. 


Time 


4/47 
Vs 
Base of Natural Logarithm 
Frror 
Angular displacement, ‘radians 
System Input 
System Output 
Time of Flight 


Angular Velocity, radians 
per second 


Time Constant 

Microfarads 

Time of Travel The total time from time 
of fire (event2) to time of 
torpedo burst (event 6). 


Time to Fvent x in 


Fig. l. 

Torpedo or Missile Speed The speed of the torpedo or 
missile with respect to the 
medium through which it is 
moving. The time rate of 
change of the distance over 
which it is moving. 

Correction to Torpedo A computed correction to 

Speed torpedo speed. 





Narre Definition 


Torpedo Turn Radius The mean radius cf curvature 
of the path taken by the torpedo 
(essentially circular) in 
traveling from the reach point 
(point 3) to the gyro involute 
point (point l). 


-* 


Target Half-Length Distance from the stern of the 
target to the middle of the target 
(MOT). 











ey oUPMERY 


An analysis of individual loops was conducted and a representative 

loop was selected and presented herein. The individual loop as 

designed will overate satisfactorily; however, using optimum settings 

of input and feedback resistances determined in section 5.1 the resrvonse 
can be improved 37.5 per cent over the best response obtained by trial 
and error, using settling time as a criteria. A further imcrovement 

was obtained by eliminating the intefral-damping network; however, 
considerable noise was incountered which would not show up in the output 
Sut would cause gear wear in the chain near the servo motor. By chang- 
ing capacitors C-3 and Cel in figure 1) to .25 microfarads instead 

of the presently installed 1 nicrofarad, the same resnonse could be 
obtained without the undesirable verformarnce at steady state. The over- 
all imorovemert in response using the resistances determined and the 


change in the capacitance was 5l;.2 per cent. 


The overall system was considered but because of the non-linear nature 

of the resolver in the circuit exact response cculd not be predicted in 
the time allowed for analysis. Linearizing assumptions were made and a 
stability analysis was performed. It was found that the system was high- 
ly stable with a gain margin of 2h db and a phase margin of 72 degrees. 
During the analysis it was discovered that the system was type zero 

which has an infinite steady state error to a ramp inrut. Within the 
area of expected performance of the system this error did not appear to 
be significant; however, if large ramp inputs are incountered or long 
c.mputing times the error can become significant. The system is easily 


compensated to becone tvpe one by an additional amplifier-motor loop. 





the system as designed it appears that by 


1 scaling that the error can be reduced significantly without 


ther compensation, 














The Gyro Lead Angle Computing Loo 


2 gyro lead anrle corputing circuit functions as a part of the 
tack Director, Mk 75. It is corsitered advantageous to familiarice 
—— 

the reader with the overall svstem prior to discussing the snecific 


circuit analvzed. 


2 Attack Directcr, Mk 75 performs as a tcrpedo data comeuter and ~isite 


ion keeper. Cortiru®:s inouts of ranre and bearing information ar surr-— 


iiied from various selected data rcwers. 
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a position keeper the Attack Jirector renerates range and Dearint af 
positi’n for transmission to sonar v:a the Analyzer Console, “» * 
the most recent values of tarret course, Sveed, rarpe, and bear ne 
Y penerates range and ~earing of apparent tarret nosition for direct 
trarsmission to sonar, fron the most recent valves cf automatically r 


wanually inserted inruts of tarret course and speed. 


As an anrle solver the Attack Director computes pvro arrle, eratiing hi, 
1 tommeco r'jnnine ‘erth fer electrically set tormedoer. 
Pequired information “or the Gertrol of wire oiided torpeices it onopl 1. 


earmec* ions to torpedo es urse ard ran @erst, and a cot. ine ore 


«< 


West? ois accowWlishéd MW Ble A) recser. 
Bio far lays of eat shiviven® 'erree ahip ere prientatel « @ ber op 
Miing the 'earin: lire to a terpet. orredd inouts are Ylgrlieed Gr 
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The Attack Director, Mk 75 functions in both automatic and display “odes. 


The specific loop which is the subject of this analvsis is a multi-loop, 
electro-mechanical servo-mechanism. The loop is comnosed of eight injivide- 
ual loops, which are interdependent in the operation of the Attack Director, 
Mk 75. Listed below are the element loops of the comvensation of lead 


angle loop and their phivsical location in the director: 


1. Gvro Angle - located in the Angle Solver Section 
2. Enabling Run = located in the Anple Solver Section 


3. Torpedo Lead Angle - located in the Position Keeper Display 
Section * 


li. Relative Target Bearing - located in the Position Keeper 
Dispvlay Section of the Attack Director 


S. Relative Tube Rearing - located in the Ballistic Computer 
Section 


6. Weapon Time to Purst - located in the Ballistic Computer 
Secticn of the Attack Director 


7. Target Speed - located in the Position Keeper Display 
Section 


8. Tarret Angle - located in the Position Keeper Displav 
Section 
The commutation of the lead angle loop is a sequential computation loop. 
This loop has four possible mechanical inputs and five possible electrical 
inputs. The outputs from this loop are three dial readings. These 
readings are valid onlv when the loop has reached a steady state condit- 


ion. 


In order to explain the operation of the computation of lead angle loop 


several assumptions were made. Initial conditions are such that Gyro 


Angle {Gj , Torpedo Lead Angle [Bo bj, Enabling Run Ic (Hm) | » and 












Weavon Time to Furst m4, are initiallv at some particular value. Another 
condition is that electric rower to all loops is not applied until such 


time as the computation begins. 


Another assumption made is that all electrical inputs are constants with 
the exception of the following variables. 

ie xm Sin Bbo 

ee Cos ets (Dmht -K1 . T26) 

Pee ym Cos i:b6 

Mee oin Bts (Dent . KI . 126) 


With the above in mind let us consider the solution to the problem: 


The firing of a preset torpedo reouires solving for a predetermined point 
of impact between tornedo and the target, noted as point 6. The torpedo 
inrut requirements are: run tc the enabling point, c(Em); the pyro anrle 
order, G3; run depth, Hvm; and a 26 Volt mark signal when O°< 1<180° fcr 


automatic horizontal search. 


Computations are made to resolve into mutually perpendicular comronents 

the final vector (leg) of the torpedo run to the tarret. This is the 
straight line connecting event l; with event 6. The comronents are ZY, 
which is parallel to the line of sight (LOS) and 3X, which is perpendicular 


to it, (reference figure 1). 


The ouantities 2 X and ZY are predicted by two isolation amplifiers summing 
the applicable inputs. These two amplifiers are the start of a seauential 
computation. The amplifier inputs involve the parallax quantities of the 
tube, Pdo and Pdn, the reach of the torpedo Rh3, the torpedo turn radius 
Ym', the relative Target Bearing B, corrected for Relative Tube Bearing 


Be, the Torpedo Lead Angle Bb6, the Target Angle Bts, target speed times 


5 








_————— 





the total run time of the torpedo (DMht . T26), and a linear spread Hé7, 
quantity if present. These quantities are related to ZX and ZY by the 


following equations. 


X = sin B (Pdo+ Rh3 cos Br - Ym' sin Br) 
- cos B (Pdn + Rh3 sin Bg = Ym' cos Bg) 


+Ym!' cos Bb6 + sin Bts (DMht . K, , (126. +67.) 


~<! 
Nt 


Rh - cos 3(Pdo + Rh3 cos Be = Ym' cos Bg) 
- sin B (Pdn + Rh3 sin Bg = Ym' cos Bg) 


- Ym' sin Bbé - cos Bts (DMht . K, . T 26 + H67) 


+ 


The cuantities £X and ZY are inputs to the two stator windings of a resol- 
ver. A servo loop positions the mechanical input of the resolver until 
one of the rotor outputs approaches null. The output of the other rotor 


is the vector sum of EX and £Y Hm6. The shaft position of the servo loop 


~4/ 2x 


is an anple equal to the tan Ey 


= Bb6, the torpedo lead arple (refer- 


ence figure 1). 


A perturbation, which includes a compensation for drift and gyro uncag- 
ing time, is summed with the null signal e(Bbé) from the resolver in an 
isolation av~plifier. The generation of this perturbation will be explained 
later. The output of this isolation amplifier is sent to a variable gain 
amplifier, VGA, which divides this signal by the other output of the lead 


angle resolver. 


The output of this amplifier is the error in Bb6. 
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This error is computed by the relation 9 = where @ is the error in Bb6, 
s is an error in the difference in length of a particular component of 
2 X and ZY and the perturbation mentioned above, and r is the vector 


sum of 2 X and = Y, Hm6. 


The vector sum of 2X and ZY is the length of torpedo travel from point 
to point 6. This signal is sent another isolation amplificr, whose 
function will be explained later. The shaft position of torpedo lead 
angle (2b6) is displayed on the outer dial of the "Target Group" and "Own 


Ship Group" dial display in the Position Keeper Display Section. 


The Bb6 shaft position is added to ~he target angle (Bts) by a mechanical 
differential (Bb6 + Bts = I). The output of this differential positions 
a limit cam, which operates a micro switch. The normally closed position 
of the microswitch gives a rirht turn mark signal to the torpedo for 
automatic horizontal search, NC, 0<1I<180 . The normally open position 
provides no mark signal and therefore, aleft turn is in order. This turn 
signal is available as a voreset torpedo function that is activated at 


the enabling point. 


The Bb6 loop also positions a two speed control transmitter synchro sy- 
stem. These CX's send an electrical analog signal proportional to the 
Torpedo Lead Angle, Bb6, to a two apeed control differential transmitter 
svnchro system that is mechaically referenced with the Relative Target 
Bearing, B. These CDX's add the two functions and yield Relative Torpedo 
Course, Cmol;. The CDX's then send this signal to the Own Ship's course 
module. The Cmols signal also references a two apeed control differential 
transmitter synchro system that is mechanically referenced with Relative 


Tube Bearing, Bg. The output of the relative tube bearing CDX's reference 











a two speed CDX svstem on the Angular Spread module sq (G). This ~odule 
is located on the firing panel. The output of these CDX's reference a 
two speed control transformer svstem on the Gyro Angle Order module, G. 


These CT's are the feedback elenent for the loop. The equation involved 


here is: IG = Cmo4— By if S4(e)] ‘ 


Bach of the terms on the right side of the above equation for G have been 
added by the use of two speed svnchro systems (CDX). This signal refer- 
ences the gyro angle module feedback CT's. These CT's have a shaft position 
propcrtional to Gyro Angle, G, and an output, which through a servo anp- 


lifier and motor, positions the gvro angle order shaft. 


The pyro angle order rear train rositions the wiper on a DC potentiometer 
(rot). This G pot is referenced with 7f7times torpedo right turn radius, 
TT Ym, and 7 times torpedo left turn radius, 7 ¥m, . The center 
tap of the pot is grounded. The voltare from the vot is positive for both 
right and left turns, but it is not of the sare magnitude, The voltage 
output at the wiper position of this pot is then proportional to the wiper 
vosition, a function of G, and the reference voltage. The analog signal 
is then the product of these two inruts; hence the arc distance traveled 


about the pgvro anple, Ym' .G. This signal is sent to an isolation 


La 
480° 
amlifier, figure 2. Other inputs to this amplifier are: The distance 

to the reach point, Rh3, from the ballistic constants network; the vector 
sum, [x X sin Bb6 + X= Y cos Bbé = Hm] , from the lead angle resolver 


previously discussed; anda correction for torpedo run depth Hvm, and 


own ship depth Hvo, Hymep 3? from the ballistic constants network. 


The output of this amplifier is proportional to Hm26. This voltage is 


sent to an impedance matching transformer in the input network of a servo 
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amplifier. The output of this network is the output of the feedback 
vernistat Hm26 on the Enabling Run servo module. The amplifier drives 
the servo motor until the two inputs are equal in magnitude and 180 out 
of phase. The motor has then driven out a mechanical analog position 
proportional to Hm26. This shaft position may be added to or subtracted 
from manually by a hand crank input labeled Fnabling Run Offset, sq(Hm). 
This is accomplished by the use of a mechanical differential. The amount 


of offset is recorded on a counter for display. 


The "nabling Run module visually displays computed Enabling Run, c(Hm), 
which is the summation of Hm26 and the offset sq(Hm). One vernistat on 
this loop sends a voltage proportional to the computed Enabling Run c(Hm), 
to the Run to Burst module as a reference source for a vernistat feedback 
loop. This particular total run to burst, [Hm26 w/e - NR = c(Hm) w/g + 
q(Hm) | » is used only in the wire puided solution and will not be discussed 


here, 


A separate cam for each weapon on the computed Enabling Run shaft, c(Hm) 
figure 3, sends a within limit signal to the gyro angle limit cam. If 
the gyro limit cams are within limits, then they send a signal to the 
depth stratum limit cams. The run depth cam, then sends a signal to the 


within limit circuitry in the weapon monitor unit. 


The Enabling Run module position six one speed control transformers that 
send an analog error signal, ec(Hm), to a servo motor in the torpedo. 
There is a control transmitter, mechanically linked to a servo motor in 


each torpedo. This CX feeds a follow-up signal oc(Hm) back to the CT's 


in the Enabling Run module to comolete the servo loop. 
The second vernistat in the Enabling Run module is referenced with the 
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difference between the sine of firing latitude and the sine of proofing 
latitude times half the product of the enabling run time and tarret turning 
rate, this minus the signal for fyro uncaging time, times own ship turn- 
ing rate. These two ocuantities are multivlied by enabling run by the use 


moetne vernistat. 
bs ot 
4 
j 


he a a f : f : 
3 é 
The output of this vernistat is sent as the perturbation to the lead an;le 


servo loop, firure ), previously discussed. 


The analog voltage, c(Hm), from the Enabling Run Vernistat, is sent to 
the total run to burst module, fifure 5, and to an isolation amp/ifier, 
which sums computed Enabling Run, c(Hm), and a torpedo path lenrth correct- 


a20n, Hmoorr 


2 
Hm = - | Sp( Hm) + Hvme (Em) + Hvmp (4m) | 
corr, Cc Cc 


(Note: The Hvme® (5) term is used only with the Mk 37 Mod O, high speed 
torpedoes.) The output of this am lifier is total torpedo run, (-) Hm27. 
This signal is sent to an impedance matching transformer in the incut 
network of the isolation amplifier of the Weapon Time to Burst servo system. 
A vernistat on the Weapon Time to Burst module is the other inrut to the 
amplifier. This vernistat is referenced with torpedo or missile speed, Um', 
from an isolation amplifier having 2? ballistic constants f(Hmv), Unt, ji(Um) 
for inouts. The f(Hvm) signal, a function of depth, is a correction to the 
torpedo or missile speed, Um. This compensation is generated on the Run 
Depth module by a pot positioned proportionally to Run Depth, Hvm, and 
referenced from the ballistic constants network with a function f(Hvm). 

The voltage output of the T26 computing vernistat is proportional to 


torpedo speed Um' times weapon time to burst T26. When the product of 
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of speed and time (ecual to total run distance Hmt), is compared with 

( - )H27 and passed through an isolation amplifier, it is an error sipnal, 
e(T26), used to position the Weapon Time to Burst, T26, servo system 
which displavs T26. A vernistat referenced with the difference between 
the sine of proofing latitude and the sine of firing latitude times half 
the torpedo turning rate and vositioned by the T26 servo system has an 
analog voltage output equal to the rroduct of weapon time to burst and the 
con-ensation for firing latitude | . T26 (sin Lyp - sin wv2)] «has 
signal is then summed at the input of an isolation amplifier with a con- 
tribution from the ballistic constant network. (-DCo . Td), where DCo 

is the Own Ship Turnine Rate and Td is the Gvro Uncaping Time. The out- 
out of the isolation amplifier is used to reference the vernistat on the 


Enabling Run module previously discussed. 


Another vernistat in the Weapon Time to Burst module sends a signal 


proportional to T26 to reference a vernistat on the Target Speed module. 


This vernistat's output is a multiple of Tarret Speed, DMht; Weapon Time 
to Burst, T26 and Kj, to yield the distance the target travels from time 
of fire to time of impact. This signal is sent to an isolation amplifier, 
figure 2, that sums this distance and linear spread quantity H67, if 
present. The linear spread inout is located in the firing nanel. Linear 
spread is never used concurrently with Angular Spread quantities, sq(G), 
previously discussed. The output of this amplifier, the total distance 
travelled, is [ (re - DMht . K,) + H67] or | (726 - DMht . Kj) + H67 | : 
and references a resolver in the Target Angle servo system. The output 
of this resolver is the sine and cosine of Tarret Angle, Bts, times the 


reference signal. These two signals become part of the inputs to the 
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(2) 


(3) 


two isolation amolifiers which vield 2X and ZY. 


Generation of one of the other inputs to each of these amplifiers bepins 
with the output of two isolation amplifiers. The inrut to these amplifi- 
ers is dependent upon the value of the Gyro Angle, C. Three conditions 
can be considered: 


(1) For 170°>G 310° and 350>G >190% the inputs into the turn- 


radius-isolation amplifier is just Ym' and the invut into the 
reach-point-isolation amplifier is Rh}3 + HVME ph 3 where Hvmgp 3 
is the correction to Rh3 for depth difference, (Hvo - Ph - Hvm). 
Under these conditions £-X and ZY are comouted as: 

ZX= | Pao sin B - Pdn cos B + (Rh3 + Hvmgpy 3) sin (B - Bre) -ym! 


cos (B = Br) + Ym' cos Bb6 + (DMht . T26 . K, + H67) sin Bts] 


For pEGre Gc 10. the incut to the pvro-turr-radius-isolation 


amplifier is zero; however, the reach-noint-isolation amplifier 


‘gs now Ph3 + Hvmgpn3 + (Ym! . G. —Z— ) . In addition 


/80° 
to this change, additional circuitry removes (-Ym" sin Bbé) 


and (Ym' cos Bb6) from the EX end ZY isolation amvlifiers. 


Hence ©=X and ZY are ccm>uted as: 


EX = Pdo sin B - Pdncos B + (Rh3 + Ym'.G. -&. ) 
“90° 
= Y = Rh - Pdo cos B + Pdn sin B + (Rh3 + Hymepy, 3 + Ym! G, 22) 


80° 
cos (B - Bg) + (Dvht . T26 . Ky + H67) cos Bts 


For 170° < G < 190° the innut to the pvro-turn-radius-isolation 
amplifier is 2 Ym' or twice the value used in case 1 above. The 


reach-point-isolation amlifier is the same as in case 1 and as 


yg 











oefore, additional circuitry removes (Ym' sin Bb6) and Ym' cos 


Bb6) from the £X and ZY isolation amplifiers. 


Furthervore, it should be noted that the quantity Ym' used in 
(1) is actually Ym} or Ym,, depending on whether the condition 


360>G >180 exists or not. 


In case (2) the same thing is true, however, its determination 


is due to the position of a not positioned in the Gyro Angle 


module, figure 6. 


In case (2) Ym' implies 2¥m'y or 2Ym'y, depending on how the Cyro 


Angle develops. 


These signals described above for case one, two, and three, reverence 
the stators on a resolver on the Relative Tube Bearing module. This 
module has two positions. One vosition is for the bearing of the 
port tube and other position is for the bearing of the starboard tude. 
The inout signals come from the own ship section where the.e two 
bearings are stabilized. The distance to the Peach Point (point 2) 
and gyro-turn radius are measured along and associated with the tube 
bearing module references these ouantities to own ship's centerline, 


(Own Ship's Course, Co.). 


To these quantities are added the appropriate parallax components, 
Pdo, Pdn, of the ship by summing the applicable input in two isolation 
amplifiers. The outnut of these two amplifiers then reference a re- 
solver on the Relative Target Bearing module. This resolver rotates 


the input quantities to distances along and across the bearing 
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line (LOS). These quantities are then sent to the summing amplifiers 


mentioned before and vield ZX and ZY. 


Objectives of Analysis 
The analvsis was undertaken to determine the response characteristics 
of the individual loops, the overall loop, check the system for stab- 


jlity, recommend potentiometer settings for optimum response and if 


found necessary recommend compensation networks. 


Period of Time Allowed for Analysis 
Work on this thesis was accomplished during eirht weeks of the summer 
period between the second and third year of the Ordance Engineering 


Guided Missiles Curriculum ending in December of 1960. 











III, Review or Introduction to Some Recent Developements in Servo 
Analysis Techniques 


3.1 Determinantal Method of Analysis, 

The theory of simple single loop control systems is well developed, and 
the block diagrams of such systems are not complicated. On the other 
hand, the block diagrams and theory of multi-loop, multi-input, milti- 
output, and multi-coupled control svstems is not very well developed 
because of the complexity of the problem involved. The block diagrams 
of such systems are complicated and difficult to intepret. One approach 
to the solution of such a problem is to formulate and solve a set of 
simultanious equations for the otek This can best be accomplished 


by arranging such systens into a standard block diagram. Arranging the 


system in this manner preatly facilitates the writing of the equations. 


Any system which has summing points in the main transmission path can 
be arranged in a block diagram by suitable block diagram manipulation 
as shown in figure 7. Writing the equations at the various nodes: 
At node a 

a= Ra - (Guy- Gagda - (G+ Gai)b - (Gog. Ggade 
A more convenient notation is: 

2 Ne a (G.,Ja - (G,)b - (Gog)c 
At node b 


Pewu, Y., & Generalized Theory of Linear @bltl-Loop wutodh tic 
Control Systems, Doctoral Dissertation, MIT, 1953 
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t node d 


d=(- Gay) d + Ra +(Gagla + Gea) d +(Gy/je 
Rearranging to bring out the symnetry 


Ut Gaal 2 ‘iGem)b 3 GC eapec + (Gda/)d = Ea 
(— Gas) 2 +lj+Grelo+ (Geode + Gaold = Pe 
PiGeae)a - Gelb eGGee)¢c tide) d + F< 


(Gag) @ - Soaid- GeaJe + UtGaasd = kd 
Jow, solving by Cramer's Rule (2 ate 


v FATAL S 
(#6 ) Gbha * Gea Gda 


- Gab C16 wo) Geb Gdb 


a = 
-Gac -Gec ec) Gde 


SEED) 8 ~.C te -Gea ({ *eaa) mas D186 we 
eee ee 
SEEM FORWAR LY AATA>S 


For simplicity let us assume the inputs 
oa = CO (2) . Ab = kai = Kd =9 


Then the output will be QO = (Gees) d 


0 b : d 
© Q@\+Ga) Gea Gea K 
B | = Gidp U+Gaa) Geb Oo 
€}- Gec Ch. “Gia@e-) “Ghee 
D= Gee tl.- Gas -Gba_—- Gea I #Gaa 
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p= Sie ("has Ki 
vA 


It is noted that AJ is equal to the characteristic equation of the sys- 
+ 


tem. For an example, take the four node system in figure 8. 


Qtdj 
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Now putting the transfer functions in the form of numerator divided by 
nominator. 
= = C2 Ly. we 
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Which when simplified becomes: 


ve es Vs Me Key 


SS 


yo <2 Ly is 4% Ky WH My es wa Wy WM yg A/ 37 


he roots of the transfer function are the zeros of the characteristic 
equation. 
bz 23 Sy os AK, Na Na Ny Ds TF Wa W3 Nv Ver =O 


ow if we desire to determine the effect of the feedback path G5, we can 


rearrange the characteristic equation in the form 


bs [bz bs De + kK, Ww M,N] 4 Wr N3 Ny WS 29 


The quantity in the brackets can be equated to zero and by root locus 
niques the roots can be determined. The characteristic equation will 

. 7} = 
then be JQ; [ (sta j(Strj Gtrs)/ tev)y +NzN3NyNMs = 9 
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The final root locus equation will then be: 


™%* M5 Ny Ns = ose (SrZ2a) es eC . be = gue 





The transfer function will then be: 


Fe ee iSe 5 
(S*#P.J (St#Has - + ee 





3.2 Determination of the Composite Open Loop Transfer Function and 


eloc ty Constant for Multi-e op ystems 
| eee I 
eo = 
l+ Fo, Yc 


| a ey a 


F< 
ST 





a = 


In terms of numerator and denominator of the closed loop transfer function 


A = Vy, 
bo - WM. 


Now applying the classical definition of the velocity constant 


Ky = Le = Ge 


SsS—~ o 


f= “eS “We 


sS~— oO Co - Mc 


It is convenient to determine the open loop transfer function in this 


manner and then apply root locus techniques. 
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_ After the open loop function has been determined, it is quiet easy to 
determine the type and order of the system. It will have the form: 
pa AAS PPO I IST 25) 
S7( 5 PAWSHALI (SAA) - ~~] 
The denominator of the open loop function is unfactored; the s" can be 
factored out prior to performing the root locus or K versus real roots 
thus the type and number of the system can be determined. If the desired 
type is not obtained, this is a criteria for selection of suitable comp- 
ensation. In any event it avoids the necessity of finding the roots of 


the system to determine the type number. 


It is often necessary to perform a root locus or K versus real root plot 
to determine the poles of the open loop transfer function. As an example, 


consider the example of section 3.1. The closed loop transfer function 
i i 
Dt bz Du Os F Ns 3 Me MS 4K, Os We Ws Ne 
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The denominator then can be factored. 
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the roots are then determined and the open loop transfer then becomes 


Sah 
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fhe closed loop roots are then determined in the usual manner. 


; Root Locus Gain versus Real Rootss 

The usual method of determining the roots of a closed loop function is to 
apply the familiar root locus principals. A somewhat similiar technique 
which embodies some of the same principals is the method of plotting 

K (root locus gain) versus the real roots. 


For illustration, let us consider this unity feedback loop shown below. 








Ot) 
> 
[The closed loop transfer function is 
(Sg, OR fo= a 
lt Fe > r Me 
The roots of Fc are where OQ, tM 20 or Ma. = Es” 
] Lb 
fhe open loop transfer function has this form: ' 
Wee 5 ade” (54 2)) AG 2a) S59 Sy) 
Do S* wth) OU tPay Os] > -- 
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2. Wheeler, R. C. H., Unpublished Class Notes 


| 


29 











se 











; plotting the poles and zeros on the S - plane: 





Fig. 9 COMPUTATION OF ROOT LOCUS | 
GAIN FROM S PLANE 


y the Evan's method the closed loop roots are determined by summing angles 

> 180° and plotting the locus of the roots. Then knowing the gain associated 
ith the system, a cut and try process is performed to locate the roots,by the 
eationship k= wet - This process is necessary along the real axis where 
00t locus exists. if these values along the axis are plotted the following 
urves will result. All the curves shown do not result from every open loop 


bransfer function but are shown as an example of all possible configeration. 
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Fig. 10 PLOT OF ROOT LOCUS GAIN 
VERSUS REAL ROOTS 









It can be shown that the maximum between /. and / is the "preakaway point" 
and the minimum between Z,and *,is the entrance point in the Evan's Root 


cus Method. 


if the system gain is given or can be chosen, a horizontal line drawn which 
intersects these curves yields the real roots of the closed loop function. 
If it does not intersect as between 4. and *% » complex roots exist. Also 
between Z, and “42 , one real root and two complex roots exist with the 


gain chosen as shown. 


si 
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This method works equally well with complex poles or zeros in the open loop 





as 


function by using the same relationship for gain K= ae < » but one 
must remember that the distances to these complex quantities to the point 


at which the gain is being computed must be included. 


the foregoing discussion, it is seen that all the real roots of the 


closed loop can be determined without plotting the root locus. 


It is a fact that the denominator of the closed loop function which is the 
characteristic equation of the system has the following form and possesses 


= 


the below listed qualities: 


ee thx) | 6FOK™ 7A --- DP =o 


= 1 

B = Sum of the roots of the closed loov function 

C = Product of the roots of the closed loop function 

Knowing the real roots from the plot of "K versus Real Roots" and using the 


relationship t.+ 22 t%3 ----- (o< ty Med FH =f me = 6 
mts +R, —->-- + Qga= 


we can now solve for the real part of the complex pair. Knowing this we can 


use the second relationship 


C, fla Sy C2 ry wt) CS gro) 2 OU 
which becomes 


We /C2z ls (2 ewe} = D 


pence «1, Wa 2 3 ANV O are known, we can solve for We thus all roots are 


known. 
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is method provides a criteria for selection of system gain since the effect 


f root location can be readily seen, and is faster than the rather 





. 


borious angle summation process necessary in the Evan's Method since in 





most cases only three points must be computed to determine a given curve. 
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IV. Derivation of Transfer Functions of Commonents and Minimun 

Velocity Constant Peauirenent 
hel Servo Arpli fier 
The frequency resvonse of the amlifier was measured and found to be flat 
to + 1 db to 7,COO0 kc. Since this is beyond the system bandwidth the 
simplification of constant gain will be used throufhout this thesis. 


It was found that this gain was avproximately 15C,000. 


Because the input signal and the various feedbacks are not summed at a 
common node (see Fig. 11), the pain could not be considered in the 
conventional manner. Rather the pain seen by the various sirnals is 
different. It was found more convenient to consider the pain as a 
function of the variable resistors. This can best be seen dy considering 
the equivalent circuit of Fig. 12. 
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Fig. 12 SQUIVALENT CIRCUIT OF INPUT NETWORK TO SERVO AMPLIFIER 
The gains as a function of the resistors Ro, Re and Ro is shown in 5.1, 
Analysis of the Individual Loop. The entire closed loop equivalent 
circuit is shown in figure 18 and sinplified block diagram equivalent 


is shown in figure 21. 
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Integral Network 


The integral-damping network is a means of compensating the amplifier 
and is shown schematically in Fig. 13. The equivalent circuit is shown 


im Fig. LW). 





Fig. 1) EQUIVALENT CIRCUIT OF INTEGRAL-DAMP NETWORK 


+ 


From Fig. li it can be shown that the transfer function is: 





Qo - S&C) 
Oo. _ (ze, C, Kw. C.)'S® + (7,.C, # &. Ca 74, GS r é 
KC, Y, aan ce 4,€C, AC 
ee Oo _ (ae) a 
Simol' fving :-—— 
Oz Kil, thilrt4#CaNS + 
~ ese rw aa 8. 


From Fig. 13 the various values of resistance and capacitance were 


obtained and substituted into the above transfer function which veilds 


a 125 
when simplified: en S®* #4255 PIO 
, ao 2. 40S 
Paes OO FSF 
Ok (s*#*sj(sr VF fo Fre We 


A frequency response analysis shown in Fig. 15 revealed that the actual 


transfer function was Go : £6.68 . This transfer function 
O« S42)(S 4200 ' 


was used in the analysis. 


It was found later in the analysis that a value of C3 = C. = .25 yf 


produced an improved response of the individual loop. Proceeding as 


IMZ3S 
13 445) (S tf //2/) 


above the resulting transfer function is . It is included 


in this section for convenience. 
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1.3 Servo Motor 
The Servo motors Mk 12, Mod O used were manufactured by KEARFOTT COMPANY, 
INC., LITTLE FALLS, NEW JERSEY. To obtain the dynamic constants, it was 


recommended that the following relationships be used. 
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B Az 
Mel VELOCITY COWSTANT)= ~ KhOrAWS PER O#O 
vt Sie lat gig ae 
VAATS 
| rc RAMANWS Ph SE*~ 
- WNeTeKk FRICTION) = Z3 7 Vi, OX E — CM. —— 
3 KAOINAS Kb SES 
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Ising a plot of number of gear teeth versus inertia in Cz#-/w* the 
total reflected inertia was found and added to the given inertia of the 


notor. The total inertia was found to be J = 3-- X/0 °° De-sn* OL 


q : ~ ee 
wy? o- 78 Gm- Cm. 


en = / Ozr-sv* 
a ee Pm 


ip TOO =~ wa IO? 
y suitable change of units this becomes: 
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When this transfer function was used in the computation to determine the 
closed loop respcnse and checked against measured response, the results 

were found to be 35 per cent in error. Several different pole-zero con- 
figerations obtained by using various values of Seelam ver feedback were 
used and the 35 per cent error was found to be consistant. It became 

obvious that the motor time constant was in error. With this in view, an 
investigation to determine the response or a second order system was attempt- 
ed by removing the tachameter and integral networks. Transient resnonse 
curves were obtained at various amplifier pain settings and found to be 


almost identical for any settinp of gain. The slight change was accredited 





to the pole of the amplifier at 6000. From these curves it can be seen 
that ts (settling time) = .355 using the relationshipt:% it was found that 
l= 44¥./4 or Ze228.2 . This can be seen readily by looking at 


the root locus plot of the second order system shown in figure 16. 
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Fig. 16 ROOT LOCUS OF INDIVIDUAL LOOP AS 
A SECOND ORDER SYSTEM 
rom the response = ee 


ie » motor transfer function then becomes: 


-. * /90 
i s€C S#2Z8.2) 


As a further check on the accuracy, a frequency response was performed. The 


It was originally 


plot of the magnitude M in decibels is shown in Fig. 17. 


intended to plot the results on a Nichols chart, transfer these results to 


a Bodie plot to find the break point or time constant. This was found too 
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practical because *here was enough backlash present to cause the phase 





anrle to shift irrationally. Using the value of Mpw from Fig. 17, second 
order curves were utilized to deternine ay the damping factor. Knowing 


§ and the natural frequency at which it occured, it was found that 


oa = « 2h 
er a ee = /3°9% 
gL: 278 


Averaging the results obtained from the transient and frequency response, 


4t was found that the motor transfer function was: 


~ 
- 


= 7 ¥.6 
(oa S( & #26) 


This transfer function was used in all the computations and was found to 


yeild excellant results. 


é 


h.2 Derivation of Theoretical “Virimum Kv Requirement 


The velocity constant reavirement was not snecified. Since the system 
under studv is a computation mechanism it was imperitive that a realistic 
minimum value of accuracv required be determined and used as a foal in 

the analysis. Within reason it would, of course, be desireable to have 

Kv as large as possible in order to reduce steady state errcrs, but at 

the same time the system need not be designed to have greater accuracy 
than sensing devices within the fire control equipment. 

¢ 

In the submarine fire control problem, bearing rate in the early stages 


‘ 5, 
of the approach will generally ranre between 0-5 per minute (rarely 
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exceeding 2 per minute). As the submarine and tarret closes, the bearing 






rate will increase, depending upon the relative motion of the two and their 
range. At very close ranges (1000 vards or less) bearing rates may become 
as high as 60° per minute, although this is rare in typical problems en- 
Peitibered by present-day high speed submarines. 

Consider the extreme case which could conceivably be encountered by two 
very-high speed submarines, each travelling 50 knots on opposite courses, 
and passing within 100C yards of one another. Their relative speed is 


100 knots (nautical miles per hour) and their range is 0.5 navtical miles. 
V= 100 Kkwors 


Y = £4000 es. 


7 


Then: Ver 8 
V = Velocity 
r = Radius 
6 = Rate of change of ® 


e Vv . 
6 = =— a = 200 Radiansf/hour 


6. 22044" 7g 8 7m: S732 
i 3600 Jéc 1@ SEC. 4B SEC. 


GC = .0555 RAOf/yr- 


NE 


















—_ OT ea, : a srt 1 ° ° ° e 
Individual components within the computer itself have accuracies generally 
of the order of 5 minutes or less. If a steady state error to a velocity 


input of the same magnitude is snecified, then from the relationship: 


vo 
cae ~ ge 
i = yr = 3./@ DEC. = 38 


ass 4 Msn. 


This would represent nearly ideal tracking under the most severe condit- 
ions likely to be incountered in a submarine fire control problem with- 
in the next decade, and is a value toward which it would be desireable 
to work. Obviously, for the great majority of applications K,, would 

not have to be nearly as large to 6btain 5 minutes or less tracking 
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V. Mathenatical Analvsis of the System and Comparison with Observed Data 





S.1 Analvsis of Individual Loops 
Since the four individual servo loops, Gyro Angle, Enabling Run, Weapon 
Time to Burst, are essentially identical, only one loop will be considered 


in detail. This is justified by measurements made on the other loops. 


in section .1, the gain cannot be considered in the conventional manner. 


: 
For illustration let us consider the Gyro Angle loop. As previously stated 
The closed loop equivalent circuit is shown in figure 18 below: 







(s #215 4200) 
K 
29% L£0% goed — GC _ a 
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| | 


Fig. 18 CLOSED LOOP EQUIVALENT CIRCUIT OF INDIVIDUAL LOOP 


It is noted that conventional summing resistors are not necessary. This 
can best be understood by considering the equivalent circuit of the input 


and feedback vernistat in figure 19. 


>See Appendix I 
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Fig. 19 EQUIVALENT CIRCUIT OF INPUT AND FEEDBACK VERNISTAT 


Further simplifying: “ 


R Kk. Oe 


70 AMIPL/elé kh 


AWD AIOTOR 


Since the reference voltage is variable and R,, > » R the circuit can 


be arranged in this manner. 


R 


— 2K e, 
79 MIMPLIEIER 
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It is convenient to derive the transfer function of the circuit comprising 
the amplifier and integral dampine network, figure 20, prior to proceeding 


with the analysis. 


26.685 
is #2) (5 r~2z200) 


GOK 


Phe 


-—— Y—» 


Fig. 20 EQUIVALENT CIRCUIT OF AMPLIFIER AND INTEGRAL DAMPING NETWORK 


Breaking the circuit in this manner is justified since neglecting the 
shunting effect of 60K and Rg = 3 negohms results in less error than 


incurred in assumed values of components. 


The resulting block diarram is: 


2a. OSS 
| (sr Zji$ 729% 





| GOA FA | 
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Simplifying: 





GOA r/aA A 2o-CS8 § 
= Sr2zJbAr200/) 






150,000 





(Cok 4 /4HZ6-68S) 
CK FIM SAAS PLO 
: i: 
Ce Se) 27 O 
Fea * closed loop transfer function of the Amplifier. 
¥ 470 
1a 150,000 (26.69S)/k 
Se a Re tke 
CALI CS4#200) 
Fe, = £:470 (S#2/) (5 fase, a 
(S42) (S4200/ (5.1) 
a |e cee tees | 26-663 
Maer fk 
The Root Locus Equation is: 
et 
( £22, 009 (26-68) Ky (a2) 
Rot 1h 
lee (St2) (S#LZ00/ 
k a (= °° ) 26-68 (/K) £ Ko00oT 2£OCuvsS CHW 
x Ae rla 


"he clese loop transfer function >f the amplifier network will ther 
2,470 (5 PA} (6 F200) 


have the form: 7a = 
4 LO*#a) (SPS) 
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The equivalent circuit of the entire loop then becomes: 3 


ea 7 See) te ae aa a ee 4 
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» 509 S Ba = ee 
I | 
R. 
A 
R | 
, 2,470 sre) (37200) 
wee (3S fib) 5( 5728) f sie 
| 20k 
(0; - Ge) ele 
———_ ee 
The block diagram representation is: 
4 20K a (2,4 TIKIPLIMSALZI(S #200) Ds) 
20K?t ky, 5S (5# 2G) (S#HGI(S#b) 
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By block diagram manipulation: 
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This simplifies to figure 21. 
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Fig. 21 REDUCED BLOCK DIAGRAM OF INDIVIDUAL LOOP 
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The root locus equation is: 
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Simplifying: 
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The closed loop function will have the form: 


Z4ESK) SHAS £290) 





__ mek» 
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It can be seen from the above equations that varying R, and Ro varies 
the location of the zero due to the tachometer and simultaneously Ro 
varies the root locus gain. It can be seen that within limits many 


settings of Ro and Ro will give identical results. 


Since Ro and R7 enter into gain as well as zero location, the conventional 
method of considering a given pole-zero confireration and varving the 

pain to determine optimum response cannot be used. Further, if a com- 
osit open loop transfer function is obtained, it can be seen that these 


factors also enter into the comoutation of the Velocity Constant, K,. 


N 
Using the relationship F, = ane the open loop function is: 
c Cc 
FEE ee, 
Denerk. (245K«) (Sr2j(strzc 





( S428//(StE)S 4200) PULA =, 5) 


y Voce ek) 
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\ 54( 20k + Ry} aN * 


LOK 
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2g a" ~ ( TateeRr ee J) (2 
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f{t is seen that three equations involving the two qualities, response and 


relocity constant, can be easily written. 
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Since the system cannot be optimmized by varying gain, a plot of these 


three equations will provide a criteria for choosing values of resistance for 


‘optimim response and maximum velocity constant. 
From figures 22 and 23 it can be seen that a near optimm Selection of 


resistance is R12 FZI0k ana Wy = 00K , It is noted that these 


values give Kv in excess of the minimm derived in l.3. 
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With these values fixed, consider the loop in three configerations: 
1. Single loop - Unity, tachometer, and integral-damp feedback 
2. Single loop - Unity, and tachometer feedback 


3- Single loop - Unity, and integral feedback 


5.2 Single Loop - Unity, Tachometer, and Integral-damp Feedback 


Substituting Ro = 300K, Ro = 100 K , and Re = 3 Megohms equation (5.1) 


becomes: 
F.. 450,000 (S#F2Z2) (S #200) 


~=—_ 
= 


cA SS#2V(SK2Z00 #13 3Y¥5 


Equation (a2), the root locus equation of the amplifier and Integral- 


damp feedback network becomes: 


SZZ5%4% 35 
ce (CS #Z/ 5 F200; 
yon 2:4¢70 (SAFZINS 42007 
“4 S* #/S 565 # ¢oO0 
ie 2,¥%70 (S421 (S #200) 
[a (S #A/5F6) (SK. 26) 


The resulting block diagram is: 
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Equation (5.3) becomes: 





(15, Joo) (S#2) (S200) 
ye = Slist2EHs £-26)(S FS 36) F SOT (SK2MSP200/j F 15, S00 GFZ NSF L200) 


Simplifying: 


 : 135, 3200 ( IF 2) (SF20e) (5.4) 
St A SGUS £452,442 5* 45,2 ¢80,6005 + 6,/20, 000 


Equation (5.4), the final root locus or characteristic equation of the 
loop is: 


zo7 (Ste) SH209 (SK YI.F/ 
=< * S( SK2BIGS FLSZE)\S H- 26) 


The root locus plot is shown in figure 2. The roots are determined from 


K versus the real roots in figure 25. 


The real roots were determined to be S = — 1862 and S = = 2.0332 . From 
the characteristic equation, denominator of equation (5.lh), the coeffi- 
cent of the second highest term is 1931 and the constant term is 
6,120,000. Proceeding as outlined in section 3.3: 
E624 2.0332 #(& fous ee A931 
Wee = 65.95 
= 33.0 


(1962j( 2. OF32) (FIAjX MFI t jH.) = 6,/29,000 
We= 22.05 


The factored closed loop function then becomes: 


15,300 (S42) (54200) | 
© S42. 0332) (S HBOS 7534 j22.65K 5433-7 22.65) 
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The response to a unit step input is: 


15,300 S#2/ (SALON fe | 
S 


axe) - e e 
(S 2.0532 KS + 1962) (St3Z347 22.65NS P33 -F 22.CS) 


The inverse Laplace trarsform yeilds: 


~ 2.03352 ~ 182 
O(é/= Le .org €E + 1.005 


= 


- "a 
- 1585 6 Gwe lhe. 5 E- (378% 


Figure 26 shows predicted © as a function of time. Figure 27 shows 


the measured response obtained. 
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302 Single Loop - Unity and Tachometer Feedback 
Removing the integral-damp network, the block diagram is: 


[+ .0O239S 


Ot) | S500 Alt/ 
S( S428) 


The closed loop transfer function is: 


15,309 
ke = 
SCS #28) F FO7( SF ¥/-9) 
Simplified becomes: 
—_ 15, ZOO 
7 S2 #395 S # 15,300 


The root locus equation is: 


ZEO7CSFHMDD 


sala sc 5*28) 


The root locus is shown in figure 28. The Root Locus Gain versus Real 


is shown in figure 29. 
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The factored closed loop transfer function is :; 


43, 300 


fo = : . 
(SF#¥¢4/ (Sh F574) 


The response to a unit step is: 


jo} 

ise 2 (+ 

O,(s) = ——=-__.____ |s] 
The inverse Laplace transform yields: 


eye ~ 3372 


a= hy t 1/32 € /32€ | ukey 


The predicted output as a function of time is shown in figure 30. The 


measured response is shown in figure 31. 


3-3 Single Loop - Unity and Integral-damp Feedback 


With no tachometer feedback, the block diagram becomes: 
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The block diagram simplifies to: 
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The closed loop transfer function then becomes: 
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Simplified: 


1/5300 \SF2) (SFL 00) 


Fe 


Bee sey s5 + 58708 SF 5,70 
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The root locus equation is: 


15, 300 ©S#2Z4NSHF200) wy 
S(S#2ONSHISIOV(SH.26/ 


The root locus plot is shown in figure 32. The plot of K versus Real 


Roots is shown in figure 33. The factored closed loop transfer function is: 


E _ 15,700 S42) \ 54200) | 
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The response to a unit step is given by: 


155200 (52) (SF260) 


ie ae 
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The inverse Laplace transform vields: 
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The predicted output as a function of time is shown in figure 3). The 


measured response is shown in figure 35. 


5.2 Qualitative Stability and Performance Analysis of Overall System 
Section 3.1 outlines the Determinantal Method of analysis used in reducing 


the functional diagram, figure 36, to the configeration shown in figure 
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It is noted that the Torpedo Lead angle Loop is inherently non-linear 


because of the resolver which is used for a function generator and a 
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SYSTEM FUNCTIONAL DIAGRAM 








nulling point. The output from the resolver is =X Cos Bb6 =- ZY Sin Bb6é 


which tends to zero when Bb6 is computed correctly. 


Because of this non-linearity, it was assumed that the inputs = X Cos Bb6 
and = Y Sin Bb6 were constant inputs and not a function of time. The 
resolver was represented as a summer. Several of the actual constants 
were combined with the scale factors of the loop for convenience and simp- 
licity. This is justified since it does not effect the accuracy of the 
stability analysis and it is highly unlikely that quantitative response 
can be predicted with the assumptions made concerning the resolver unless 


a prohibitive number of small perturbations are considered. 


The assumption is also made and justified in section 6.1 that for signals 
which do not saturate the amplifier or motor that the knabling Run, 
Weapon Time to Burst, and Gyro Lead Angle transfer functions are essentia- 
lly identical. For large signals this assumption is especially invalid 

in the Weapon Time to Burst Loop because of the 40 to 1 gear ratio which 


permits the motor to saturate. 


With these assumptions in mind, the "Standard Block Diapram" becomes 


that shown in figure 37. Where: 
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It is convenient to expand the delta so that the ro can be obtained 
since it will be used in the numerator of the closed loop function. 
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The closed loop function becomes 
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Substituting the various values for the numerator and denominator of the 
functions, the open and closed loop transfer functions become: 
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Observing equation 5.5 the open loop transfer function is of the form 


Fo» A(S#2) (SAH) where n = QO. This indicates, of course, that the 
SO  SMSHP IC SAHA) 
system is type zero. 
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The steady state error to a ramp is: 
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This indicates that the steady state error would be infinite and the sy- 


stem must be compensated in such a way to make it a type one if ramp 


inputs of large magnitude or long duration are expected to be incountered. 
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Since, of necessity, certain linerizations were made; therefore, no 
conclusions can sensibly be made concerning stability except that the 
necessary condition that all the coefficients of the characteristic 
equation are positive is met. This being the circumstance, a frequency 
response analysis was made. The results were plotted on a Nichol's 
chart. It was found that there was a gain margin of 26 db and a phase 
margin of fe This, of course, shows that the overall system is highly 
stable. The results obtained from the Nichol's chart were then plotted 
on a Bodie diagram, figure 38, to determine the dominant time constant 

of the system and to acertain experimentally that the system was actually 


type zero. These results are discussed in section 5.3-b. 


5.3 Discussion of Analysis ami Observed Data 


5.3a Individual Loop 
The individual loop was analyzed in three possible configerations,unity, 


tachometer, and integral-damp feedback; unity and tachometer feedback; 


and unity and integral-damp feedback. . 


By using figure 22 and 23 it was found that near optimum values of 
resistance were R, = 300 K and Ro = 100 K. It was found that the value 
of Re used had little effect on the response of the system; however, 
considerations of noise in the system dictated a value of approximately 
three megohms. Using these values of resistance, computations were made 
to predict the response of the loop. The predicted and measured response 


is shown in figure 26 and 27 respectively. 


Prior to the determination of the values of resistances, the best response 


obtained was a critically damped case with a settling time in excess of 
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cl) seconds. It can be seen that the predicted response and measured 
response agree exactly and a settling time of ~15 seconds was obtained. 


This shows an improvement in response of 37.5 per cent. 


The analysis of the loop without the integral-damp network was performed 
in section 3.2. The measured response is shown in figure 31. Comparing 
this with the predicted response in figure 30, it is again seen that the 
measured and predicted values agree. It is also noted that an improvement 
in response was obtained with a settling time of .11 seconds. This,of 
course, is desireable; however, at steady state there was audible noise 
caused by gear chatter. This noise did not show up in the output but 
would cause excessive wear in the gear train near the motor; therefore, 
the integral-damp network is considered desirable. It was found analy- 
tically that by changing capacitors C3 and C), (See figures 13 or 1h) to 
25M a predicted and measured response with a settling time of .12 
seconds was obtained without any associated noise problems. The theo- 


retical pole location is calculated in section l).2. 


The last configeration considered involved unity and integral-damp feed- 
back. From figures 3) and 35 it can be seen that the predicted and 


measured response are again in agreement except the peak overshoot is 





somewhat greater in the actual system. Settling time is .3 seconds. 
It is noted that substantial tachometer feedback must be used since for 


smooth computation each loop must be critically damped. 


5.3b Overall System 
As discussed in section 5.2, the overall system is highly stable with a 


phase margin of i and a gain margin of 2h db. The data from the Nichol's 


plot is used to construct figure 38. From this it is noted that the 
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first asymptote has a zero db slope; therefore, type zero system. The 
time constant agrees reasonably with the actual obtained from the meas- 
ured response shown in figure 39. The theoretical settling time is .6 
seconds and the measured is .7 seconds (See figure 39). The response 


to a ramp input is shown in figure 10. 


Since the theoretical computations and frequency response indicated that 
the system was type zero with a finite error to a step input and an 

infinite error to a ramp input, it was considered advantageous to meas- 
ure these errors. The inputs and resulting errors are shown in figures 


hl and 2. : 


The results are therefore conclusive that the system is in fact type zero. 
While the mlti-loop system is unmodified (i.e., a type O servomechanism), 
there will be a constant error to a step and a corresponding, increasing 
error to a ramp input. The measurements in figures }1 and 2 show the 
magnitude of the errors involved. The first is a measurement of the error 
versus time for the step and the second is the measurement of the error 
versus time corresponding to the ramp. From the first the error ratio 


to the step input can be calculated to be: 





+ OV 2 . coeey 
540475 


i 
In a typical submarine problem, the total bearing change from initial 


contact to the firing point is approximately 150 degrees. This means 


that, under these circumstances, the error at firing point will be: 


(.0053#) (50°) = . 8° 


One degree = 17.6 yards/1000 yards. Therefore, the following errors 
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FIGURE h1 


MEASURED ERROR OF THE OVERALL SYSTEM TO A STEP INPUT 
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FIGURE 2 


MEASURED ERROR OF THE OVERALL SYSTEM TO A RAMP INPUT 





world result with a computed error in the Gyro Lead Angle Computation 
Loop of 0.8": 
at: @,000 yds. : 1.05 yas. 

2,000 yds. : 28.10 yvds. 

3,000 yds. : 2.15 yds. 

l,,000 yaSe : 56.20 yds. 
These errors are, of course, in addition to any other errors generated 
in other parts of the computer. It can be readily seen that if fire 
control problems of short range with the resultant large ramp or long 
ranges where tracking is done for long periods of time, a significant error 


may develope. 


~ 


Looking at equation (5.5),the open loop transfer function of the system, 
it is seen that an interration is required or an "S" is needed in the 
denominator to have a type one system. It appears that a feed forward 


loop added to the block diagram, figure 37, would accomplish the desired 


results. The resultant standard bleck diagram is shown in figure 13. 





Fig. 3 STANDARD BLOCK DIAGRAM WITH COMPENSATION 
NECESSARY TO PRODUCE A TYPE ONE SYSTEM 
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The determinant for delta becomes 
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The closed loop function is 
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It can be seen that the necessary and sufficient condition for a type one 


system has been met since oa; zs SCS*2) can be factored from each tern. 
@ 
Where Go ® a is of the form: Gg = AT 
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VI.e Conclusions and Recommendations 


6.1 Single Loop 
The single loop designed and constructed by Librascope Division of General 


Precision Corporation will produce satisfactory response without further 
changes in design. However, by judicious selection of the various input 
resistors, an improvement is possible of 37.5 percent using settling 
time as a criteria. It can be seen from the root locus, figure 2), and 
K versus Real roots, figure 25 that the system is absolutely stable for 
all values of gain. It can also be seen that the gain can vary appreci- 
ably without effecting the response. Should the gain become very low, 
an appreciable residue will be associated with the root near the zero 

at S = <2. Tachometer feedback is essential to provide a critically 
damped system for computation. Integral-damp feedback slows down the 
response by 26.5 per cent from that obtained without it; however, the 
damping is required to prevent high frequency oscillations in the gear 
train during steady state operation. In selecting values of capacitance 
(C, - C), = .25 4#F ) in figure 1, the determental effect on response 
due to the integral-damp feedback is eliminated. The recommended changes 
improves the response a total of 5.2 per cent over the best obtained 


by trial and error on the physical system. 


6.2 Overall System 


The overall system is type zero as designed. This is satisfactory for 
most applications intended since the error involved is small; however, 
should large ramp inputs or tracking problems of long duration be encounte- 
ered, it is possible to generate excessively large errors. The solution 


to this problem is a modification involving a motor and amplifier in a 
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feed forward loop. This loon will be physically connected at the error 
input to the Torpedo Lead Angle and feed forward to the output of Bbé. 


This is represented in the modified standard blcck diarram, figure l2. 


An alternate method which may be prefered, since the primary function 

of the loop is to mechanize an equation, is to determine a mean ramp, 
calculate the various outpouts as a function of time, and scale the varicus 
loops to reduce or eliminate the resultant error. A check must, of course, 
be made over the full range of ramps expected, after the scaling has been 
accomplished, to insure that the error does not become excessive else- 


- 
Zz 


where. 


Measured frequency and transient resoonse indicates that the system operation 
will be satisfactory with .7 seconds settling time to a unit step, gain 


margin of 2) db, and phase margin of 7. 
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VII. Recommendations for Further Investigation 


Should it be deemed desirable to predict accurately the response to a 
complex multi-loop system utilizing resolvers, it would be advantageous 

to use the linearizations made herein, consider small pertubations, and 
determine their accuracies obtained bv measurements on the physical system. 
This would provide a criteria for evaluating future analysis which may be 


conducted during design work on other systems. 


In addition it may be advantageous to consider linearization technicues 


other than the ones used herein. 


If it is required that the overall svstem be optimized, an analogue computer 


study is warranted. 
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APPENDIX I 
Induction Resolver 


An induction resolver is an electromagnetic device employed as 

a Dasic element in analog computers, automatic control systems, 
angle data transmission systems and plan position indicating 
radars. It is used for performing trigonometric computations, 

The induction resolver is essentially a transformer with rotary 
variable coupling between primary and secondary windings (Figure 1). 
The windings are distributed in the various stator and rotor slots 
in such manner as to obtain a true sinusoidal relationship between 
output voltage and angle of rotation of the rotor. Since the 
units are used at audio and power frequencies, it is customary 

to use a high-permeability laminated core similar to that used in 


an induction motor. 


— wn > 





guts SIV 9 
Fig. 1 - Basic Resolver 
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It is customary to treat the stator winding as the primary or 
input winding, and the rotor winding as the secondary or output 
winding. If, as in Figure 2, there are two secondary windings 
at right angles to each other (in "space quadrature"), one sec- 
ondary will have induced in it a voltage which is proportional 
to the sine of the angle of rotation between it and the primary, 
and the other secondary will have induced in it a voltage which 
is proportional to the cosine of the angle of rotation between 


it and the primary. 
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Fig. 2 = Resolver with two Rotor Windings 
To make induction resclvers even more useful, two windings are 
provided on both the stator and rotor. The pairs of stator wind- 
ings and rotor windings are located at right angles to each 
other as shown in Figure 3. This permits computation of cosines 
as well as sines, and even more complex problems involving both 
sines and cosines can be solved when both stators are excited. 
In addition to the resolution of input voltages into sine and 
cosine components, resolvers can perform the synthesis of resul- 


tants from components and the rotation of coordinate systems. 
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Fig. 3 - Resolver Input and Output Relationship 


As shown in Figure lL, currents in the two stator windings 
set up magnetic fields at right-angles to each otrer to form a 
resultant flux vector. This flux vector may assume any angular 
position, depending upon the amplitudes and polarities of the 
input currents. Its magnitude and proportional to the vector 
sum of the instantaneous values of the input currents, and its 
angular position is equal to the arctangent of the ratio of 
the two currents, For AC input voltages within the operating 
frequency range of the resolver, the rotor voltages are propor- 
‘_.tional to the magnitude of the flux vector and the sine or cosine 
of the angle of rotation of the rotor winding with respect to 
the vector. The output voltages may be expressed in terms of 
the input voltages and the rotor angular position, by the fol- 


lowing equations: 
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The subscripts refer to the various rotor and stator windings 


and ©& is the rotor angular position. 





ed 


>< Fig. h - Stator Flux Vectors, 
While other types of operations are possible with a resolver, 
its use has been principally with three types of inputs, namely; 
(1) sine wave voltages of the same electrical phase, (2) sine 
wave voltages shifted 90 degrees apart, and (3) irregular wave 
shapes, principally of fips sawtooth or square wave variety. For 
ie first type of operation, the output voltages remain fixed 
in electrical phase shift and their magnitudes vary in accordance 
with equations (1) and (2). For the second and third types of 
operation, the output voltages may also be determined by equations 
(1) and (2) if the instantaneous values of the voltages are taken 
into consideration. For the case of the resolver with two inputs 
90 degrees apart in electrical phase, each of the two rotor output 
voltages remains constant in amplitude, independent of rotor an- 
gular position. However the phase of the voltages varies contin- 
uously with rotor angular position, The phase shift between the 


two output voltages remains fixed at 90 degrees. In this mode of 
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operation, then, the equations for the output voltages are: 
&e= & Ll. (3) 


¢ 


fy = & Le t¥o" (1) 


For the third case, where irregular wave shapes are used, op- 
eration most frequently falls into two categories. In phase, 
square waves or flat topped sine waves are used in some computing 
applications. The performance is adequately described by equations 
(1) and (2) without modification if the two input wave shapes 


are the same, which is usually the case. 


Sawtooth wave shapes are frequently used in radar sweep appli- 
cations. In the simplest case, a sweep voltage is applied to 
one stator and the two output voltages from the rotor windings 
are the same sweep voltage but with an amplitude proportional 
to the sine or cosine of the rotor angular position. If these 
two voltages are applied to the "X" ard "Y" deflection plates, 
or coils of a cathode ray tube, they cause a radial trace to be 
formed. The trace will assume the same angular position as the 
resolver rotor. This type of operation finds application in 


plan position indicator types of radar displays. 
While the arrangement shown in Figure 5a is suitable for schema- 


tic representation of a resolver, the arrangement shown in Figure 


Sb is more suited for block diagrams and functional schematics. 
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The diagram, equivalent schematic and equations are shown for 
clockwise rotation (facing the shaft end) and apply to all 
standard types of resolvers. Vvhere it is desired to change the 
direction of rotation of the resolver, it is only necessary to 
interchange input A with B and output C with D. This can be 


noted on the block diagram by a negative sign precedinge . 
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A resolver used in the simplest form of computation is show in 
Figures 6a and b. One stator winding is excited with a voltage 
proportional to the hypotenuse of a right triangle and the rotor 
is positioned to the angle @ . The two output voltages are pro- 
portional to the two legs of the triangle. 
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Fig. 6 - Elementary Resolver Application 
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In Figure 7 the problem is similar, but the hypotenuse and one 
of the legs, Y are known. A servo amplifier and motor are used 
tq position the resolver rotor until the voltage, ER os is del- 
ivered by the resolver equals the known voltage Y. The voltage 
PR, is then proportional to the length of the other leg of the 
triangle, and the angle © is the angle between the second leg 


and the hyrotenuse. A solution of this type may be used in de- 


termining elevation angle and ground range when slant range and - 
depth are known. x= fe, , Y7ER, 
REF. y-/~ Pur R= fs, , OF Sim Y 
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Fig. 7 - Angle Computation 


In Figure 8, the inputs are two vectors 90 degrees apart. The 
servo positions the resolver rotor until one output is zero. 

At this position the other rotor being 90 degrees away, has maxi- 
mum coupling with the flux vector formed by the X and Y imputs. 
Its output is therefore proportional to the vector sum of these 
two inputs. The rotor angular position is the angle between one 


leg and the hypotenuse of the triangle. 





Fig. 8 = Vector Addition 
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In the rotation of coordinates, both stators and rotors are 

used, With rectangular coordinates ¥ and Y as inputs to a resol- 
ver, as in Figure 9, the outputs are also rectangular coordinates 
of the same point referenced to a new set of axes rotated with 
respect to the first. The operation if fully described by the 
equations in Figure 3 which will be recognized as those required 


for coordinate rotation. 


The actual operation can best be understood by breaking it into 
two steps. The X-Y coordinates impressed on the two stators 

form a flux vector in the resolver. This vector is independent 

of rotor position. Its angular position is the arctangent of Y/X. 
Each of the rotor windings develops a voltage proportional to 

the sine or cosine of its angular position with respect to the 
flux vector (not with respect to zero). If the reference is made 
with respect to zero then the outputs follow the equations 

in Figure 9. “he rotor voltages are the X', Y' components of 

the flux vector. The angle between the X', Y' coordinate ref- : 


erence of the inputs is the same as the rotor anguler position. 


my, 
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Fig. 9 - Coordinate Rotation 
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INDUCTION POTENTIOMETER 


An induction potentiometer, or linear synchro transmitter, may 

be considered a special kind of resolver which provides 

accurate linear indication of shaft rotation about a reference 
position in the form of a polarized voltage whose magnitude is 
proportional to angular displacement, and whose phase relation- 
ship indicates direction of shaft rotation. The principal differ- 
ence between a resolver and an induction potentiometer is that 

the latter's output voltare varies directly as an angle and not as 
the sine of that angle. Its rotation is usually limited to less 


than 290".. 





Fig. 1 - Induction Potentiometer Voltage Output 
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These devices are analogous to resistance potentiometers, but 

since they are induction type components they have less re- 
straining forces acting upon their rotors, and hence are capable 
of providing better resolution. Naturally, induction poten- 
tiometers do not require sliders such as those used in resistance 
tipes. Therefore, circuit interruptions are eliminated, no wear 
occurs as a result of rubbing parts, and accuracy is consequently 
continously maintained at the original level throughout the 
operational life of these types of components. These advan- 
tageous features are virtually 2 necessity in certain applications. 
For example, in gyroscope systems where low restraining torques 
and low pick-off angular errors are required, wide use is made 

of these devices. In general, induction potentiometers find use 
in applications where resistive potentiometers are impractical, 
principally because of the following features: (1) induction 
potentiometers, having no wiping contacts, may be used as 
gyroscope pick-offs since they contribute less spurious friction 
torque: (2) input and output are isolated: (3) resolution is 
infinite: (lh) noise level is low: and (5) the total angle of travel 


is limited to less than 180° 


Constructed very much like a resolver, induction potentio- 
meters differ in that their windings and slots are not uniformly 
distributed, but instead are deliberately modified to produce a 


linear output. Normally they have only one excitation winding 
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and one output winding, the former usually carried on the rotor 


and connected by means of two sliprings and two brushes. 


The use of induction potentiometers as linear pick-offs may 
be considered as a special resolver application in which mutual 
coupling varies not as Z,,Siw@, but as ZO. With zero input im- 


pedance and infinite output impedance, the voltage gradient, 


£ 
H, may be defined as H ‘i - Under this condition the trans- 


fer function may be expressed as: 


fe = ee. AO ee f~ ma F 
L 7% Z>5 i o* Zp 


where 7s = secondary impecance 


Zz, = primary impedance 


NN 
" 


. = load impedance 


Z.s input impedance 
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